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We predict the features of the Collins function, which describes the fragmentation of a transversely 
polarized quark into an unpolarized hadron, by modeling the fragmentation into pions at a low 
energy scale. We use the chiral invariant approach of Manohar and Georgi, where constituent 
quarks and Goldstone bosons are considered as effective degrees of freedom in the non-perturbative 
regime of QCD. To test the approach we calculate the unpolarized fragmentation function and the 
transverse momentum distribution of a produced hadron, both of which are described reasonably 
well. In the case of semi-inclusive deep-inelastic scattering, our estimate of the Collins function in 
connection with the transversity distribution gives rise to a transverse single spin asymmetry of the 
order of 10%, supporting the idea of measuring the transversity distribution of the nucleon in this 
way. In the case of e + e~ annihilation into two hadrons, our model predicts a Collins azimuthal 
asymmetry of about 5%. 



I. INTRODUCTION 

The influence of transverse spin and transverse mo- 
mentum on fragmentation processes is at present a 
largely unexplored subject. The Collins fragmentation 
function Jl), correlating the transverse spin of the frag- 
menting quark to the transverse momentum of the pro- 
duced hadron, could give us the first chance to study this 
effect. Moreover, being chiral-odd, the Collins function 
can be connected to the transversity distribution func- 
tion, which is chiral-odd as well, and thus can allow the 
measurement of this otherwise elusive property of the nu- 
cleon, which carries valuable information about the dy- 
namics of confined quarks. Beside being chiral-odd, the 
Collins function is also time-reversal odd (T-odd). 

In spite of the apparent difficulty in modeling T-odd 
effects, in a recent paper Q we have shown that a non- 
vanishing Collins function can be obtained through a con- 
sistent one-loop calculation, in a description where mas- 
sive constituent quarks and pions are the only effective 
degrees of freedom and interact via a simple pseudoscalar 
coupling. 

In our previous work [^j little care has been devoted to 
the phenomenology of the Collins function. In contrast, 
our interest here lies in obtaining a reasonable estimate 
of this function and the observable effects induced by it. 
At present, only one attempt to theoretically estimate the 
Collins function for pions exists (U , and little phenomeno- 
logical information is available from experiments. The 
HERMES collaboration reported the first measurements 



of single spin asymmetries in semi-inclusive DIS [^,^| , giv- 
ing an indication of a possibly non-zero Collins function. 
The Collins function has also been invoked to explain 
large azimuthal asymmetries in pp 1 ' — > irX |BQ . In this 
case, however, the extraction of the function is plagued 
by large uncertainties, due to the possible presence of 
hadronic effects both in the initial and final state, and 
hence does not allow any conclusive statement yet. Re- 
cently, a phenomenological estimate of the Collins func- 
tion has been proposed combining results from the 
DELPHI, SMC and HERMES experiments. However, in 
spite of all the efforts to pin down the Collins function, 
the knowledge we have at present is still unsufficient. 

In this work we calculate the Collins function for pions 
in a chiral invariant approach at a low energy scale. We 
use the model of Manohar and Georgi |9| , which incorpo- 
rates chiral symmetry and its spontaneous breaking, two 
important aspects of QCD at low energies. The spon- 
taneous breaking of chiral symmetry leads to the exis- 
tence of (almost massless) Goldstone bosons, which are 
included as effective degrees of freedom in the model. 
Quarks appear as further degrees of freedom as well. 
However, in contrast to the current quarks of the QCD 
Lagrangian, the model uses massive constituent quarks - 
a concept which has been proven very succesful in many 
phenomenological models at hadronic scales. With the 
exception of Ref. ||Tcf| , the implications of a chiral in- 
variant interaction for fragmentation functions into Gold- 
stone bosons at low energy scales remains essentially un- 
explored. To investigate the Collins function for vector 
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mesons like the p lyj is beyond the reach of the approach. 

Although the applicability of the Manohar-Georgi 
model is restricted to energies below the scale of chi- 
ral symmetry breaking A x w 1 GeV, this is sufficient to 
calculate soft objects. In this kinematical regime, the 
chiral power counting allows setting up a consistent per- 
turbation theory [ fL^ ] . The relevant expansion parameter 
is given by l/A x , where I is a generic external momen- 
tum of a particle participating in the fragmentation. To 
guarantee the convergence of the perturbation theory, we 
restrict the maximum virtuality p 2 of the decaying quark 
to a soft value. We mostly consider the case p? = 1 GeV 2 . 

The outline of the paper is as follows: We first give the 
details of our model and present the analytical results of 
our calculation. Next, we discuss our results and com- 
pare them with known observables, indicating the choice 
of the parameters of our model. Then, we present the 
features of our prediction for the Collins function and its 
moments. Finally, using the outcome of our model, we 
estimate the leading order asymmetries containing the 
Collins function in semi-inclusive DIS and in e + e~ anni- 
hilation into two hadrons. 



II. CALCULATION OF THE COLLINS 
FUNCTION 

Considering the fragmentation process of a quark into 
a pion, q*(k) — > n(p)X, we use the expressions of the 
unpolarized fragmentation function D\ and the Collins 
function in terms of light-cone correlators, depend- 
ing on the longitudinal momentum fraction z of the pion 
and the transverse momentum k T of the quark. The def- 
initions read llSjl^ 



D 1 {z,z 2 kl)= j- z J dfc+Tr[A(fc,p) 7 -] 
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with mjr denoting the pion mass and = e w (we 
specify the plus and minus lightconc components of a 
generic 4- vector a p according to a ± = (a ± a 3 )/-\/2)- 
The correlation function A(k,p) in Eqs. ([j], ^), omitting 
gauge links, takes the form 



*Note that this definition of slightly differs from the 
original one given by Collins Q. 
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We now use the chiral invariant model of Manohar and 
Georgi |J to calculate the matrix elements in the corre- 
lation function. Neglecting the part that describes free 
Goldstone bosons, the Lagrangian of the model reads 



C = 4> + Y - m + 9aAib) 1> ■ 



(4) 



In Eq. (g) the pion field enters through the vector and 
axial vector combinations 
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with u = exp(i r-7r/2F w ), where the Tj are the generators 
of the SU(2) flavour group and = 93 MeV represents 
the pion decay constant. In absence of resonances, the 
pion decay constant determines the scale of chiral sym- 
metry breaking via A x = 4atF w . The quark mass m and 
the axial coupling constant g A are free parameters of the 
model that are not constrained by chiral symmetry. The 
values of these parameters will be specified in Sec. III. 



Though we limit ourselves here to the SU(2) flavour 
sector of the model, the extension to strange quarks is 
straightforward, allowing in particular the calculation of 
kaon fragmentation functions. For convenience we write 
down explicitly those terms of the interaction part of the 
Lagrangian that are relevant for our calculation. To 
be specific we need both the interaction of a single pion 
with a quark and the two-pion contact interaction, which 
can be easily obtained by expanding the non-linear rep- 
resentation u in terms of the pion field: 
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Performing the numerical calculation of the Collins func- 
tion, it turns out that the contact interaction (Q), which 
is a direct consequence of chiral symmetry, plays a dom- 
inant role. 

At tree level, the fragmentation of a quark is modeled 
through the process q* — + irq, where Fig. |l| represents the 
corresponding unitarity diagram. Using the Lagrangian 
in Eq. (^|), the correlation function at lowest order reads 
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This correlation function allows to compute the unpo- 
larized fragmentation function D\ by means of Eq. (Q), 
leading to 
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FIG. 1. Lowest-order unitarity diagram describing the 
fragmentation of a quark into a pion. 



(b) 



Note that the expression in (||) is only weakly dependent 
on the transverse momentum of the quark. In fact, D\ is 
constant as a function of k T , if = and (or) m = 0. 
Because our approach is limited to a soft regime, we will 
impose an upper cutoff on the k T integration, as will be 
discussed in more detail in Sec. [II. This in turn leads to 



+ h.c. 



a finite D\(z) after integration over the transverse mo- 
mentum. 

The SU(2) flavour structure of our approach implies 
the relations 

= D{-^ = Df^ w ~ = Df**~ = 2 D x , (11) 
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FIG. 2. One-loop corrections to the fragmentation of a 
quark into a pion contributing to the Collins function. The 
hermitian conjugate diagrams (h.c.) are not shown explicitly. 
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FIG. 3. One-loop self-energy, and vertex corrections. 



where D\ is the result given in Eq. (|9|). In the case of 
unfavoured fragmentation processes D\ vanishes at tree 
level, but will be non-zero as soon as one-loop corrections 
are included. According to the chiral power counting, 
one-loop contributions to Di are suppressed by a factor 
I 2 /A 2 compared to the tree level result. The maximum 
momentum up to which the chiral perturbation expan- 
sion converges numerically can only be determined by an 
explicit calculation of the one- loop corrections. 

Like in the case of a pseudoscalar quark-pion cou- 
pling j|], the Collins function turns out to be zero 
in Born approximation. To obtain a non-zero result, we 
have to resort to the one-loop level. In Fig. || the corre- 
sponding diagrams are shown, where we have displayed 
only those graphs that contribute to the Collins function. 
The explicit calculation of if -j 1 is similar to our previous 
work g]. The relevant ingredients of the calculation are 
the self-energy and the vertex correction diagrams. These 
ingredients are sketched in Fig. || and can be expressed 
analytically as 
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where flavour factors have been suppressed. For later 



purpose, we give here the most general parametrization 
of the functions S, T± and T 2l 



E(jfc) = Afc + Bm, 



(15) 



IMM = j§rlz\C l +D l 1>+E l $ + F l l>^ , (16) 

T 2 (k,p) = ^-75 (C 2 + D 2 i> + E 2 Jt + F 2 i> ft) . (17) 

The real parts of the functions A, B, C\, D\ etc. could 
be UV-divergent and require in principle a proper renor- 
malization. Here, we don't need to deal with the question 
of renormalization at all, since only the imaginary parts 
of the loop diagrams are important when calculating the 
Collins function B. 

Taking now flavour factors properly into account, the 
contributions to the correlation function generated by the 
diagrams (a) , (b) and (c) in Fig. ^ are given by 
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The correlation functions of the hermitian conjugate dia- 
grams follow from the hermiticity condition A^?'(fc, p) — 

7 °A; i) (fc,p)7 . 

Summing the contributions of all diagrams and insert- 
ing the resulting correlation function in Eq. (Q) , we even- 
tually obtain the result 
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Thus, the Collins function is entirely given by the imagi- 
nary parts of the coefficients defined in Eqs. Q15- 17). We 



can compute these imaginary parts by applying Cutkosky 
rules to the self-energy and vertex diagrams of Fig. |^. 
Explicit calculation leads to 
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where we have introduced the so-called Kallen function, 
A(fc 2 , m 2 , m 2 ) = [k 2 — (m + m 7r ) 2 ][fc 2 — (m — m^) 2 ], and 
the factors 

h= ( d 4 l6(l 2 ~ml)6((k-l) 2 -m 2 ) 
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These integrals are finite and vanish below the thresh- 
old of quark-pion production, where the self-energy and 
vertex diagrams do not possess an imaginary part. 



Thus, Eq. ^ in combination with Eqs. (j2|)-((2^) 
gives the explicit result for the Collins function in the 
Manohar-Georgi model to lowest possible order. Because 
of its chiral-odd nature, the Collins function would van- 
ish in this model if we set the mass of the quark to 
zero. The same phenomenon has been observed in the 
calculation of a chiral-odd twist-3 fragmentation func- 
tion [^0). The result in Eq. ( pl| ) corresponds, e.g., to 
the fragmentation u — > 7r°. The expressions of the re- 
maining favoured transitions are obtained in analogy to 
Eqs. ( pO| , pT[ ). Unfavoured fragmentation processes in the 
case of the Collins function appear only at the two-loop 
level. 



III. ESTIMATES AND PHENOMENOLOGY 

A. Unpolarized fragmentation function and the 
choice of parameters 

We now present our numerical estimates, where all re- 
sults for the fragmentation functions in this subsection 
refer to the transition u — > tt + . To begin with we calcu- 
late the unpolarized fragmentation function Di(z) which 
is given by 



Di(z) 



Ki 



dK 2 D^z.Kl), 



(27) 



where K T = —zk T denotes the transverse momentum of 
the outgoing hadron with respect to the quark direction. 
The upper limit on the K 2 integration is set by the cutoff 
on the fragmenting quark virtuality, /i 2 , and corresponds 
to 
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= z (1 — z) /J 2 — z m 2 — (1 — z) m 2 . 



(28) 



Besides m and g A , the cutoff /i 2 is the third parameter 
of our approach that is not fixed a priori. However, as 
will be explained below, the possible values of p? can be 
restricted when comparing our results to experimental 
data. Unless otherwise specified, we always use the values 



0.3 GeV, 
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1 GeV 
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At the relevant places, the dependence of our results on 
possible variations of these parameters will be discussed. 
Few remarks concerning the choice in ( p9] ) are in order. 
The value of m is a typical mass of a constituent quark. 
The choice for the axial coupling can be seen as a kind of 
average number of what has been proposed in the litera- 
ture. For instance, in a simple SU(6) spin-flavour model 
for the proton one finds g A « 0.75 in order to obtain the 
correct value for the axial charge of the nuclcon |J. On 
the other side, large iV c arguments favour a value of the 
order of one |Tq ], while, according to a recent calculation 
in a relativistic point-form approach jig], a g A slightly 
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above one seems to be required for describing the axial 
charge of the nucleon. Finally, our choice for fi 2 ensures 
that the momenta of the outgoing pion and quark, in the 
rest-frame of the fragmenting quark, remain below val- 
ues of the order 0.5 GeV. In this region we believe chiral 
perturbation theory to be applicable, meaning that our 
leading order result should provide a reliable estimate. 

In Fig. |^ we show the result for the unpolarized frag- 
mentation function D\^^ . Notice that in general the 
fragmentation functions vanish outside the kinematical 
limits, which in our model are given by 




corresponding to the situation when the upper limit of 
the K 2 integration becomes equal to zero. We con- 
sider our tree level result as a pure valence-type part 
of . The sea-type (unfavoured) transition u — > 

7r + is strictly zero at leading order. Therefore, we 
compare the model result to the valence-type quantity 
j-)u— >ir _ ; where the fragmentation functions 

have been taken from the parametrization of Kretzer frif | 
at a scale Q 2 = 1 GeV 2 |] Obviously, the ^-dependence of 
both curves is in nice agreement. We point out that such 
an agreement is non-trivial. For example, in the pseu- 
doscalar model that we used in our previous work 0, D\ 
behaves quite differently and peaks at an intermediate 
z- value. 

On the other side, we underestimate the parametriza- 
tion of Ref. by about a factor of two. Some remarks 
are in order at this point. Although a part of the discrep- 
ancy might be attributed to the uncertainty in the value 
of g A , the most important point is to address the ques- 
tion as to what extent we can compare our estimate with 
existing parametrizations. The parametrization of Jl7| ] 
serves basically as input function of the pQCD evolution 
equations, used to describe high-energy e + e~ data, and 
displays the typical logarithmic dependence on the scale 
Q 2 . A value of Q 2 — 1 GeV 2 is believed to be already 
beyond the limit of applicability of pQCD calculations. 
On the other side, our approach displays, to a first ap- 
proximation, a linear dependence on the cutoff /i 2 . It is 
supposed to be valid at low scales and it is also stretched 
to the limit of its applicability for fi 2 = 1 GeV 2 . In 
this context it should also be investigated to what ex- 
tent the inclusion of one-loop corrections, which allow 



T Other parametrizations ]18|-|2Cj use a starting energy scale 
Q 2 > 2 GeV 2 , which is too high to allow a comparison with 
our results. 
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FIG. 4. Model result for the unpolarized quark fragmen- 
tation function D"^" (solid line) and comparison with a 
parametrization of Ref. [[L7| (grey line). 




FIG. 5. Model result for the average hadron transverse mo- 
mentum for different choices of the cutoff: /i 2 = 0.5 GeV 2 
(dotted line), ^ 2 = 1 GeV 2 (solid line), ^ 2 = 1.5 GeV 2 
(dashed line) and comparison with a fit to experimental re- 
sults from DELPHI jn) (grey line). 

for the additional decay channel q* — > 7T7r<7, will increase 
the result for D x at p 2 = 1 GeV 2 . Finally, we want to re- 
mark that to our knowledge there exists no strict one-to- 
one correspondence between the quark virtuality /i 2 and 
the scale used in the evolution equation of fragmentation 
functions, which in semi-inclusive DIS, e.g., is typically 
identified with the photon virtuality Q 2 . For all these 
reasons, a smooth matching of our calculation and the 
parametrization of |l7| cannot necessarily be expected. 
Despite these caveats, the correct z-behaviour displayed 
by our result for D\ suggests that the calculation can 
well be used as an input for evolution equations at a low 
scale. In the next subsection we will elaborate more on 
this point in connection with the Collins function. 

The best indication of the appropriate value of the cut- 
off /x 2 may be obtained when comparing our calculation 
to experimental data of the average transverse momen- 
tum of the outgoing hadron with respect to the quark, 
which we evaluate according to 

k 2 

(\K T \)(z) = [ AKl^D^Kl). (31) 

Uiyz) Jo 

In Fig. [3 we show the result of this observable as a func- 
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tion of z for three different choices of the parameter p 2 . 
As a comparison, we also show a fit (taken from Ref. gj) 
to experimental data obtained by the DELPHI collabo- 
ration ^T|. Like in the case of -Di(z), the shape of our 
result is very similar to the experimental one, which we 
consider as an encouraging result. For p 2 = 1 GeV 2 our 
curve is about 30% below the data. Such a disagreement 
is not surprising, keeping in mind that at LEP energies 
higher order pQCD effects (e.g. gluon bremsstrahlung, 
unfavoured fragmentations, etc.) play an important role, 
leading in general to a broadening of the K T distribu- 
tion. For experiments at lower energies, however, where 
pQCD contributions can be neglected in a first approxi- 
mation, it may be possible to exhaust the experimental 
value for {\K T \)(z) with genuine soft contributions as de- 
scribed in our model. This in turn would determine the 
appropriate value of the cutoff /i 2 . For example, such 
a method of matching our calculation with experimen- 
tal conditions could be applied at HERMES kinematics, 
even though the method is somewhat hampered since K T 
is not directly measured in semi-inclusive DIS. In this 
case, one rather observes the transverse momentum of 
the outgoing hadron with respect to the virtual photon, 
PhA_, which depends on both K T and the transverse mo- 
mentum of the partons inside the target p T . At leading 
order in the hard scattering cross section one can in fact 
derive the relation 

2 tt j &p 2 T p 2 T h{x,p 2 T ) ttJ dKlKlD^z^Kl) 
fl(x) D^z) 
= z 2 (p 2 T )(x) + (K 2 )(z) , (32) 

where x represents the Bjorken variable. 

B. Collins function 

We now turn to the description of our model result 
for the Collins function. In Fig. [| is plotted for 
three different values of the constituent quark mass, 
m = 0.2, 0.3, 0.4 GeV. In a large z-range, the func- 
tion does not depend strongly on the precise value of 
the quark mass, if we choose it within reasonable limits. 
That's why we can confidently fix m = 0.3 GeV for our 
numerical studies. It is very interesting to observe that 
the behaviour of the unpolarized fragmentation function 
D\ is quite distinct from that of the Collins function: 
while the former is decreasing as z increases, the latter 
is growing. 

The different behaviour of the two functions becomes 
even more evident when looking at their ratio, shown in 
Fig. We emphasize that also from the experimental 
side there exists some evidence for an increasing ratio 
Hi/D\. In a recent analysis of the longitudinal sin- 
gle spin asymmetry measured at HERMES, Efremov et 
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FIG. 6. Model result for the Collins function for different 
values of the constituent quark mass: m = 0.2 GeV (dotted 
line), m = 0.3 GeV (solid line), m = 0.4 GeV (dashed line). 
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FIG. 7. Model result for H^/Dl 

al. Q extracted a behaviour H^/D\ cx z for z < 0.7. 
We consider the agreement in finding a clearly rising ra- 
tio I Di as remarkable, even though in the analysis of 
Ref. [|| some simplifying assumptions have been used in 
order to obtain information on from data taken with 
a longitudinally polarized target. It will be very inter- 
esting to see if dedicated future experiments can confirm 
such a behaviour. We also mention that ratios of the 
Collins function or any of its moments with D\ are almost 
independent of the coupling constant g A . The reason is 
that the one-loop correction containing the contact inter- 
action is only proportional to g\ , as D\ is, and is domi- 
nating on the others. Furthermore, the ratio / D\ is 
nearly independent of the cutoff p 2 . In conclusion, the 
prediction shown in Fig. |7] is almost independent of the 
choice of parameters in our approach. 

At this point we would like to add some general re- 
marks concerning the z-behaviour of our results. It turns 
out that the shape of all the results does not vary much 
when changing the parameters within reasonable limits. 
In particular, variations of g A and of the cutoff p 2 only 
change the normalization of the curves but not their 
shape. In this sense our calculation of fragmentation 
functions has a strong predictive power. This has a di- 
rect practical consequence if one uses, for instance, our 
result of the Collins function as input in an evolution 
equation: the ^-dependence of the input function can be 
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FIG. 8. Model result for H 1 L(1/2) fLh (solid line) and com- 
parison with the product ((\K T \) /2zm 7T ) (H^ / Di) (dashed 
line). Note that the positivity bound requires the ratio to be 
smaller than 0.5. 



0.4 




0.2 0.4 0.6 0.8 1 



FIG. 9. Model result for H HX/2) /Di (solid line) and com- 
parison with the same ratio, where H^ 1 ^ 2 ^ is calculated ac- 
cording to Eq. @ with Mc = 0.3 GeV (dashed line) and 
M c = 0.7 GeV (dotted line). 



adjusted to the shape of our H^ , while its normalization 
can be kept free in order to account for uncertainties in 
the values of g A and /z 2 . 
In Fig. p| we plot the ratio 
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which enters the transverse single spin asymmetry to 
be discussed in the following subsection. This quantity 
rises roughly linearly within a large z-range, leading to 
a similar z-behaviour of the transverse spin asymmetry. 

is no longer independent of the cutoff /i, 
but rather the same dependence as in the case of ( | K T | ) 
(shown in Fig. ||) can be assumed. In Fig. |[ this ratio is 
compared to the expression 
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2zm- K Dt(z) 77 Df(z) 
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A very close agreement between the two different curves 
can be observed, indicating that the model predicts 
a quite similar transverse momentum dependence of 
both the Collins function and D%. In the literature, 
this feature is sometimes assumed in phenomenological 
parametrizations of H^. Note, however, that in our ap- 
proach deviations from this simple behaviour can be ex- 
pected, if Di is also calculated consistently to the one- 
loop order. 

The Collins function has to fulfill the positivity 
bound I22U23II 



(35) 



^L H t(z,K 2 T )<l Dl (z,K 2 T ) 



Integration over K 2 gives the simplified expression 

Cr4 < 36 > 

Di(z) 2 



which is satisfied by our model calculation. It is clear, 
however, that increasing the value of /i 2 will eventually 
result in a violation of the positivity condition. To avoid 
such a violation, we should calculate D\ and consis- 
tently at the same order, i.e., the one- loop corrections to 
D\ should be included. By doing so, the positivity bound 
will be fulfilled even at larger values of ^ , for which our 
numerical results are no longer trustworthy. 

From our results, we expect an increasing behaviour of 
the azimuthal asymmetry in p^p — + tt X as function of 
Xp, qualitatively similar to what has been predicted in 
Ref. H in the context of the Lund-fragmentation model. 
At this point, it is also interesting to discuss the com- 
parison of our results with the ones obtained using the 
so-called "Collins guess". On the basis of very general 
assumptions, Collins suggested a possible behaviour for 
the transverse spin asymmetry containing [Q. This 
suggestion has been used in the literature (see, e.g., 
Refs. [p^-p7[) to propose the following shape for the 
Collins function 



Ht (h \z) = v 



\K T \ M c 



2z 



Ml 



A'; 



■Di (37) 



with the parameter Mq ranging between 0.3 and 
0.7 GeV. Using our model outcome for the unpolarized 
fragmentation function, we apply Eq. ( |37| ) to estimate 
Hi , and in Fig. iwe show how this compares to 
the exact result of Eq. (|33|). There is a rough qualitative 
agreement with the Collins ansatz for the lowest value 
of the parameter Mc, although it is not growing fast 
enough compared to the exact evaluation. On the other 
side, in the Manohar-Georgi model there is no agreement 
with the Collins ansatz for high values of the parameter 
Mq, which might indicate that the relation suggested in 
Eq. ( j37| ) should be handled with care. 

Finally, we display in Fig. UU the quantity 



H 



-L(l) 



Di(z) D 1 {z) 



2z 2 m' 



Ht{z,K 2 T ), (38) 
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FIG. 10. Model result H^ W /D! (solid line) and compari- 
son with the product ((K%) /2z 2 m%)(Hi / ' Di) (dashed line). 



because this ratio appears in the weighted asymmetries 
to be considered below. In Fig. 10, also the expression 



(Kl){z)Ht{z) Ht(z) 

— IT 



2z 2 ml D^z) 



D{(z) 



K 2 

dK 2 a t n ^ 



T 2z 2 m1 



D 1 (z,K J T ) 



(39) 



d 6 CT[/ T + dVuj. 




T_ 

2 



Ph 



(40) 



where the subscript U indicates an unpolarized electron 
beam, the index a denotes quark flavours, and /i is 
the usual unpolarized quark distribution in the nucleon. 
Subtracting the cross sections we obtain the polarized 
part 0, 



d (Jui — d 6 aui = —\Sj 



4a; 



sxy 



2 sin(0f 



d 2 p T d 2 k T S 2 pi 



P, 



Ph± ■ k T 
\Ph±\m„ 



Y,e 2 a K(x,P 2 T )Ht a (z,z 2 k 2 T ). 



(41) 



is shown for comparison. Once again, there is a remark- 
able agreement between the two different expressions, 
confirming the quite similar K T behaviour of H^- and 
Di. 



Integration over the azimuthal angles would cause the po- 
larized part of the cross section to vanish. After defining 
the angle <j> = e/>f — 2<f>f , we consider the sine/> weighted 
transverse spin asymmetry 



(sine 



C. Asymmetries in semi-inclusive DIS and e e 
annihilation 

We turn now to estimates of possible observables con- 
taining the Collins function. We will take into consider- 
ation one-particle inclusive DIS, where the Collins func- 
tion appears in connection with the transversity distri- 
bution of the nucleon, and e + e~ annihilation into two 
hadrons belonging to two different jets. 

In the first case, we consider the DIS cross section with 
a transversely polarized target^] and the production of 
one pion. We denote the transverse polarization vector 
of the target as St- The cross section is differential in six 
variables, for which we choose x, y, z, \Ph±\, 4»h > <Pi > where 
Phi_ is the transverse component of the pion momentum, 
<frh is its azimuthal angle with respect to the target spin, 
and <pf is the azimuthal angle of the lepton scattering 
plane again with respect to the target spin. 

Orienting the spin of the target in two opposite di- 
rections and summing the cross sections we isolate the 
unpolarized part Jl4j|, 



* Transverse vectors and azimuthal angles are defined as 
lying on a plane perpendicular to the direction of the virtual 
photon. 



= J djf d 2 P h± Sine/. (dW T - d 6 er ax ) 
/ dcf>? d 2 P h± (d 6 aui + dVi) 



(42) 



Inserting Eqs. ([ll 40) into the definition of the asym- 
metry results in an expression where the transverse mo- 
menta of hi and are still entangled in a convolution 
integral ^8|. To resolve the convolution, it is required 
to assume a particular dependence of the transversity 
distribution on the intrinsic transverse momentum. The 
simplest example is 



hi{x,p 2 T ) « hi(x) 



5{p 2 T ) 



(43) 



that means supposing there is no intrinsic transverse mo- 
mentum of the partons inside the target. Under this as- 
sumption, the pion transverse momentum with respect to 
the virtual photon is entirely due to the fragmentation 
process, i.e., Ph± = Kt = —zkx, and the convolution 
can be disentangled 



(sin 4>) UT (x,y,z) 



5* (1-1/) E a elhi(x)Ht {i)a (z) 



xy \ 



(44) 



where the approximation sign reminds that the equality 
is assumption-dependent. 
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If we want to disentangle the convolution integral of 
Eq. (41) without making any assumption on the intrinsic 
transverse momentum distribution, we need to weight 
the integral with the magnitude of the pion transverse 
momentum p9[ . This procedure results in the azimuthal 
transverse spin asymmetry 



IF/ 



sin^) (x,y,z) 



UT 



J dcf>? d 2 P h± ^ sin0 ( dVc/T - d 6 a vl ) 



f 



if d 2 A±(dV c/T + d 6 ^) 

_L(1) 0/ 



\S T \ 



(45) 



We achieved an assumption-free factorization of the x 
dependent transversity distribution and the z-dependent 
Collins function. The measurement of this asymmetry 
requires to bin the cross section according to the mag- 
nitude of the pion transverse momentum. On the other 
side, this asymmetry represents potentially the cleanest 
method to measure the transversity distribution together 
with the Collins function. Moreover, it is not afflicted by 
complications due to Sudakov factors |50|. 

We show predictions for both transverse spin asymme- 
tries defined in Eqs. and ((45]). Different calculations 
can be found in the literature, e.g., in Rcfs. |U 32,27]]. To 
estimate the magnitude of the asymmetries, we need in- 
puts for the distribution functions, in particular for the 
transversity distribution. Several model calculations of 
this function are available at present (see |33j for a com- 
prehensive review). We refrain ourselves from consider- 
ing many different examples and rather restrict the anal- 
ysis to two limiting situations. In the first case we adopt 
the non-relativistic assumption h\ = gi, while in the sec- 
ond case we exhaust the upper bound on the transversity 
distribution, i.e., hi < \{fi + gi) HJ]. We use the simple 
parametrization of gi and /i suggested in [[35). At the 
moment, more sophisticated parametrizations are avail- 
able, taking also scale evolution into account. However, 
all these parametrizations are compatible with each other 
to the extent of our purpose here, that is to give an es- 
timate of the asymmetries for a low scale. We focus on 
the production of n + , where the contribution of down 
quarks is negligible, not only because of the presence of 
unfavoured fragmentation functions, but also because the 
transversity distribution for down quarks is supposed to 
be much smaller than for up quarks. 

In Fig. ^ we present the azimuthal asymmetry de- 
fined in Eq. ( [44] ) as a function of x, after integrating 
numerator and denominator over the variables y and z, 
for the two cases described above. In Fig. |l2l we present 



In Fig. _ 

the same asymmetry as a function of z, after integrating 
over y and x. As already mentioned before, our predic- 
tion is supposed to be valid at a low energy scale of about 
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FIG. 11. Azimuthal transverse spin asymmetry (sin <j))uT 
as a function of a;. Solid line: assuming hi = gi. Dashed line: 
assuming hi = + <?i). The functions fi and gi are taken 
from 01. 
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FIG. 12. Azimuthal transverse spin asymmetry (sin 4>)ut 
as a function of z. Solid line: assuming hi = gi- Dashed line: 
assu ming hi = \{fi +gi)- The functions fi and gi are taken 
from |3ql. 



lGeV 2 . Neglecting evolution effects, it could be utilized 
for comparison with experiments at a scale of few GeV 2 . 
We assume the value of the transverse polarization to 
be \St\ = 0.75. In performing the integrations, we apply 
the kinematical cuts typical of the HERMES experiment, 
as described in Therefore, our prediction is partic- 
ularly significant for HERMES, which is supposed to be 
the first experiment to measure this asymmetry. In prin- 
ciple, the simultaneous study of the x and z-dependence 
of the asymmetry yields separate information on the dis- 
tribution and fragmentation parts and allows to extract 
both up to a normalization factor Note, however, 

that this procedure relies on the assumption of up-quark 
dominance and is valid only if the asymmetry is truly 
factorized, so that the x-dependence can be ascribed en- 
tirely to the distribution functions and the z-dependence 
entirely to the fragmentation functions. Kinematical cuts 
could partially spoil this situation. We would like to 
stress that our calculation predicts an asymmetry up to 
the order of 10%, which should be within experimental 
reach, and suggests the possibility to distinguish between 
different assumptions on the transversity distribution. 
Using the same procedure as before, we have esti- 
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FIG. 13. Azimuthal spin asymmetry ( ^- L sin^)j/T as a 
function of x. Solid line: assuming h\ — g\. Dashed line: 
assu ming hi = + gi)- The functions /i and gi are taken 
from 1331. 
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FIG. 15. Azimuthal asymmetry {P%j_ cos 2(j>} E 
annihilation into two hadrons, integrated over the range 
0.2 < Z2 < 0.8 (solid line), and over the range 0.5 < Z2 < 0.8 
(dashed line). 
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FIG. 14. Azimuthal spin asymmetry ( ^- L sin (f>)uT as a 
function of z. Solid line: assuming hi — g±. Dashed line: 
assu ming hi — + <7i). The functions /i and g\ are taken 
from 13511. 



mated the asymmetry defined in Eq. (jig), containing 
the weighting with \Ph±\/Tn-n- The results are shown in 
Fig. [l3| as a function of x and in Fig. |lj as a function 
of z. The magnitude of this asymmetry is higher than 
in the unweighted case, which is due to the fact that the 
weighting spuriously enhances the asymmetry by about 
a factor two. 

Besides appearing in semi-inclusive DIS in connection 
with the transversity distribution of the nucleon, the 
Collins function can be independently extracted from an- 
other process, that is electron-positron annihilation into 
two hadrons belonging to two back-to-back jets p6| , |37| | . 
We restrict ourselves to the case of 7-exchange only. In 
this process, one of the two hadrons (say hadron 2) de- 
fines the scattering plane together with the leptons and 
determines the direction with respect to which the az- 
imuthal angles must be measured. The cross section is 
differential in five variables, e.g., z\, Z2, y, \Ph±\, 4>- The 
variables z\ and z 2 are the longitudinal fractional mo- 
menta of the two hadrons. In the center of mass frame 
y = (1 + cos0)/2, where 8 is the angle of hadron 2 with 
respect to the momentum of the incoming leptons. The 
vector Ph± denotes the transverse component of the mo- 



mentum of hadron 1 and <j) is its azimutal angle with 
respect to the scattering plane. For a more detailed de- 
scription of the kinematical variables we refer to J36],[37| . 
We define the azimuthal asymmetry 

<p h \ COS 20) e+e _ 

= / d 2 P h ±P% ± cos 20 d 5 a e+e - 

2 sin 2 9 H^ (1) ( Zl ) H^ {1 \z 2 ) 

- l + co S 2 (D lizi) D{ 1) (z 2 )+D[ 1 \z 1 )D 1 (z 2 ) 

where summations over quark flavours are understood. 
The weighting with a second power of Ph± in the numer- 
ator is necessary to obtain a deconvoluted expression. 
We prefer to use the same weighting in the denominator 
as well, to avoid a modification of the asymmetry just 
caused by the weighting. 

In Fig. [l5] we present the estimate of the asymmetry de- 
fined above, entirely based on our model. The asymmetry 
has been integrated over z 2 and 9, leaving the dependence 
on z\ alone. We have extended the 6 integration interval 
all the way to [0,7r], to obtain a conservative estimate. 
In fact, limiting the interval to [j, ^-] will enhance the 
asymmetry by a factor 2, approximately. Because the 
Collins function increases with increasing z, we also get 
a larger asymmetry by restricting the integration range 
for z%. As an illustration of this feature, in Fig. [15] we 
present two results, obtained from two different integra- 
tion ranges. Our prediction is supposed to be valid only 
at low energy scales and should be evolved for compar- 
ison with higher energy experiments. It is important to 
note that we estimate the asymmetry to be of the order 
of about 5%, and thus should be very well observable in 
experiments. 
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IV. SUMMARY AND CONCLUSIONS 

We have estimated the Collins fragmentation function 
for pions at a low energy scale by means of the Manohar- 
Georgi model. This model contains three essential fea- 
tures of non-perturbative QCD: massive quark degrees of 
freedom, chiral symmetry and its spontaneous breaking 
(with pions as Goldstone bosons). Because of the chi- 
ral invariant interaction between pions and quarks, the 
fragmentation process can be evaluated in a perturbative 
expansion. The constituent quark mass, the axial pion- 
quark coupling g A and the maximum virtuality /i 2 of the 
fragmenting quark are free parameters of our approach. 
The quark mass and g A are constrained within reasonable 
limits. To ensure the convergence of the chiral perturba- 
tion expansion, /i 2 cannot exceed a typical hadronic scale. 
We have mostly considered the value /i 2 = 1 GeV 2 , which 
guarantees that the momenta of the particles produced 
in the fragmentation process stay well below the scale of 
chiral symmetry breaking, A x w 1 GeV. To determine the 
appropriate value of /i 2 , the average transverse momen- 
tum of a data set could be used. In any case, we observed 
that variations of the free parameters within reasonable 
limits have only a minor influence on the shape of the re- 
sults, implying that our approach has a strong predictive 
power for the z-behaviour of the various functions. 

We have found that the Manohar-Georgi model repro- 
duces reasonably well the unpolarized pion fragmenta- 
tion function and the average transverse momentum of a 
produced hadron as function of z, supporting the idea of 
describing the fragmentation process by a chiral invariant 
approach. 

Compared to the unpolarized fragmentation function, 
modeling the Collins function is considerably more diffi- 
cult, mainly because of its chiral-odd and time-reversal 
odd nature. In our approach, the helicity-flip required 
to generate a chiral-odd object is caused by the mass 
of the constituent quark, while the T-odd behaviour is 
produced via one-loop corrections. The Collins function 
exhibits a quite distinct behaviour from the the unpo- 
larized fragmentation function. In particular, the ratio 
H^-/Di is strongly increasing with increasing z. 

On the basis of our results, we have calculated the 
transverse single-spin asymmetry in semi-inclusive DIS 
where the Collins function appears in combination with 
the transversity of the nucleon. This observable will be 
measured in the near future at HERMES and could also 
be investigated at COMPASS, Jlab (upgraded) and eR- 
HIC. For typical HERMES kinematics the asymmetry is 
of the order of 10%, giving support to the intention of 
extracting the nucleon transversity in this way. We be- 
lieve that our estimate of the Collins function, despite its 
uncertainties, can be very useful for this extraction. 

More information on the Collins function from the ex- 
perimental side is urgently required. In this respect, the 



most promising experiment seems to be e + e~ annihila- 
tion into two hadrons, where appears squared in an 
azimuthal cos 20 asymmetry. According to our calcula- 
tion, an asymmetry of the order of 5% can be expected, 
which should be measurable at high luminosity accelera- 
tors, such as BABAR and BELLE. Dedicated measure- 
ments of the Collins function would be extremely impor- 
tant for the extraction of the transversity distribution. 
Moreover, they could answer the question whether a chi- 
ral invariant Lagrangian can be used to model the Collins 
function. 
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